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We analyse the effect of synchronization between noise and periodic signal in a two-state 
spatially extended system analytically. Resonance features are demonstrated. To have the 
maximum cooperation between signal and noise, it is shown that noise strength at resonance 
should increase linearly with the frequency of the signal. The time scale of the process at 
resonance is also shown to increase linearly with the period of the signal. 
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Noise induced large response in bistable potential to a weak periodic signal conimonly referred to as 
stochastic resonance (SR) attracts considerable interests during recent years |]l]-^. Manifestation of SR is 
usually exihibited as non-monotonous behavior of power spectrum of the process or as a similar behavior 
of the peak heights of the escape-time distribution related to switching between two wells, as a function of 
the noise strength [Q-0. The usage of the term "resonance" is then criticised based on two facts [||. The 
value of the noise strength for which the maximum of the spectral density at forcing frequency occurs does 
not simply correspond to the period of the signal and for fixed noise strength, the spectral density decreases 
monotonically with frequency of the signal, at variance with the notion of ordinary resonance. One therefore 
would naturally be curious whether any characterisation of this phenomena in terms of our usual definition 
of resonance is possible as it was originally pictured [pi . 

Recently, it has been found that such description is indeed possible [foj. Simulating the continuous 
stochastic process into a stochastic point process with the help of an analog circuit they llfl] measure the 
strength of the first few peaks of the residence time distribution. It is found that when strength of the 
first peak is plotted as a function of the frequency of the signal keeping noise strength constant, it hits a 
maximum for that period which is equal to twice of the inverse of the Kramers rate for the unbiased process 
for that particular noise strength. Alternately, when strength of the first peak is plotted as a function of 
noise strength it attains a maximum at that value of the noise strength for which inverse of Kramers rate 
for the unbiased process is equal to the period of forcing signal. Subsequently, this resonant behavior is also 
experimentally verified in the polarised emission of vertical cavity surface emitting laser pT| . In this paper 
we analyse these results analytically. 
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FIG. 1. Potential V{x,t) at time t as a function of x. 

Response of a two-state spatially extended system embedded to a noisy environment under the infiuence 
of a periodic field is of our concern. The simplest model |l2| for the two-state spatially extended system is 
to consider a particle moving in the piecewise double well potential V{x) under the influence of white noise. 
The potential V{x) is shown in Fig.l. The influence of the external periodic field is usually described by the 
modulation of the potential well in the following fashion. The potential at any instant of time t for the left 
state is replaced by — | Vb | +Acosil.t and that for the right state is replaced by — | Vb | —Acosil.t with A, H. 



being the amplitude and frequency of the periodic signal and | Vo | is the barrier height when there is no 
modulation. 

The Fokker-Planck equation (FPE) for the probability distribution function P(x, t) for position x of the 
particle at time t for this process is 
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where D is strength of the white noise. It is clear that for the potential in Fig.l, ^ = everywhere except 
at the discontinuous points. Therefore in each region of V{x, t) ^constant, the FPE (1) reduces to simple 
diffusion equation. The solution in each region are to be matched with the continuity of probability current 
and jump conditions at the discontinuous points [— fca/2, ka/2] at each time. These conditions are to be 
supplemented by the reflecting boundary conditions at the walls [— a(l + fc/2), a(l + fc/2)]. 

In this paper we however are interested in first passage time density function (FPTDF) with mean first 
passage time (MFPT) being a important parameter of the process. To be specific we concentrate on the 
events starting from x — —1 and ending at a; = 1. As the potential is symmetric, the MFPT < t{—l ^ 1) > 
would be same as < t{l — > —1) >. Thus the conditions that the required solution of Eq.(l) satisfies are 

P'(-a(l + fc/2),t) = , (2a) 

e[''\v«\+Acosnt]/Dp^_j^^l2 - 0, t) = P{-ka/2 + 0, f) , (2b) 

P'{-ka/2 - 0, t) = P'{-ka/2 + 0, t) , (2c) 

P{ka/2 - 0, t) = e[-l^''l-^'=°''^*l/-°P(fca/2 + 0, t) , (2d) 

P'{ka/2 - 0, t) = P'{ka/2 + 0, i) , (2e) 

P(l,i) = 0, (2f) 

where prime over P in the above equations denote derivative with respect to x. Eq.(2a) corresponds to 
reflecting boundary condition for the wall at x = — a(l + k/2), Eq.(2b) and Eq.(2d) are the jump conditions 
at the discontinuous points x = ±fca/2, Eq.(2c) and Eq.(2e) are the matching condition for continuity of 
probability current and Eq.(2f) is the absorbing boundary condition ar x = 1 for the process considered. 
For (A/ D) small, we expand the solution in each region as 

oc 

P(x,t) =e[-w-A*-(A/i5M*)]^(^)[l+ ^(^/^)™/„^(^^i)] ^ (3) 

m— 1 

with /o = 1, and g{t) being defined as 

g{t) ^ cosflt ; -a(l + fc/2) <x < -ka/2 , 

= ; -fca/2 < a; < fca/2 , (4) 

= —cosflt ; fca/2 < a; < 1 . 

In Eq.(3), A is some constant and fm{x, t) are functions, which are to be determined. Substituting the ansatz 
(3) in Eq.(l) and equating equal powers of (A/D) we obtain the following set of equations in each region: 



dt dx"^ 



^^ = -{\/Dmx) , (5a) 

[ilj(x)fn{x,t)]= g{t)fr,-i{x,t) ;n = 1,2,... , (5b) 



with /o = ij^ix). For small {A/D) we keep only Eq.(5a) which is obtained after equating the terms 0(1). 

Substituting the approximate form of probability density in Eqs.(2a)-(2f) we obtain the following set of 
matching conditions for the functions tpi^)- 

V';(-a(l + fc/2)) = 0, (6a) 

e[-l^«ll/-°?A^(-fca/2 - 0) = ^j^,{-ka/2 + 0) , (6b) 

V';(-fca/2-0)=^;(-fca/2 + 0), (6c) 

V'^(fca/2 - 0) = e[-l^''ll/-°?A^(fca/2 + 0) , (6d) 

V^;, (fca/2 -0)=7A; (fca/2 + 0), (6e) 

^m(I) = , (6f) 



We note that the solutions of Eqs.(6a)-(6f) are the exact eigenfunctions of the corresponding unperturbed 
process with associated eigenvalues A^. Thus with this approximation the conditional probability for getting 
the particle at position x at time t when it is known to start from x — —1 ati = Ois 

Pix, t I -1, 0) = e-[^(-)-^(-i)l/2^e-^[s(*)-il/^ J2 e'^"*V'n(x)Vn(-l) , (7) 

n 

where the value of V{x) and g{t) would be taken as the respective values of the region where x falls. As tpn 
form the complete set, P{x,t ~ 0+ | —1,0) = S{x + 1) is automatically satisfied. The Eqs. (6a)-(6f) are 
readily solved to obtain the full set of eigenfunctions and corresponding eigenvalues. They express as 

■>pf,{x) = C^cos{kf,[x + a{l + fc/2)]} ; -a(l + fc/2) < a; < -fca/2 , (8a) 

= -Cp[el^''l/2-°sm(fc^a)sm{fc^(x + fca/2)} 

_e-l^"l/2^cos(fc^a)cos{fc^(a; + fca/2)}] ; -fca/2 < a; < fca/2 , (8b) 

~ C^[cos{k^ka)cos{k^[x + a(l — fc/2)]} 

— e' "'' sin{kfj,ka)sin{kfj,a)cos{kf^{x — fca/2)} 

—e^'^"'' ^ sin{k^ka)cos{kp,a)sin{kf_t{x — fca/2)}] 

; fca/2 < a; < 1 , (8c) 

where constant C/^ is determined from the normalisation condition 

Clj ^4'l{x)dx^l. (9) 

The corresponding eigenvalues are determined from the transcendental equation 

cos{kf^ka)cos{k^[l + a(l — fc/2)]} 
— e' "'' sin{kfj_ka)sin{kfj_a)cos{kf^{l — fca/2)} 
_e-l^«l/'f^sm(fc^fca)cos(fc^a)sm{fc^(l - fca/2)} = . (10) 

The expressions (8) and fc^ obtained from Eq.(lO) are substituted in Eq.(7) with A^ = Dk^ to obtain the 
conditional probability P{x,t \ —1,0) at any point x at time t. 

The survival probability at time t , S{t) can be obtained from Eq.(7) as 



5(0. /; 



P{x,t\-l,{))dx , (11) 



and the first passage time density function (FPTDF) g{t) is obtained from S{t) as 

.W--^. (12) 

Physically, g{t)dt gives the probability that the particle reaches x = 1 for the first time in the time interval t 
and t + dt starting from x = —1 at i = 0. As the expression of the conditional probability (7) is approximate 
and strictly valid for small (A/D), we should demand the necessary physical conditions like positivity, 
normalisability on g(t). Different moments of g(t) can be readily calculated from the normalised g{t) as 



<t^ >= 



/ t^g{t)dt ; J = 1,2,... . (13) 
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FIG. 2. FPTDF g{t) as a function of f for frequency 
n = 0.6 and (D/Vo) = 1.0, 2.0, 2.7, 3.2. 
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FIG. 3. Peak areas as a function of (D/Vo) and 
MFPT < i > as a function of (D/Vo) with signal 
{A — 0.2) [line with filled circles] and without signal 
{A — 0) [line with open circles]. 



The calculations are done with the potential parameters, a = 1, fc = 1, | Vq \~ 0.25, A — 0.2. Normalised 
FPTDF g{t) as a function of t for fixed value of signal frequency fi = 0.6 and different values of {D/ | Vq |) 
are plotted in Fig. 2. Multiple peaks are observed for small values of noise strength and as noise strength 
increases background is found dominant. g(t) extends for a longer time for small noise strength. The peak 
heights fall exponentially with time for fixed {D/ | Vq |) but peak heights are increased as noise strength is 
lowered. Successive peaks are separated by a time exactly equal to the period of the signal frequency showing 
that the probability of transition is maximum after each period. At these particular times coperation between 
noise and the signal in making a transition is more. It clearly demonstrates that the signal synchronizes the 
noise induced transition. Peak areas (counting peaks from left ) are also found to be decreasing with time. 

On subtracting the exponential background, the peak areas P„ = J g{t)dt, where the integration is done 
around n th. peak are plotted as a function of (D/ | Vo |) in Fig. 3. The figure demonstrates a non-monotonous 
behavior for each peak exihibiting a signature of maximum coperation or synchronization between noise and 
external signal at specific value of noise strength for a particular signal frequency fi. The position of the 
maximum of P„ shifts towards lower value of noise strength as peak number increases. All curves merges to 
low D showing at very low noise strength transition would mainly be controlled by signal most of the time. 
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FIG. 4. FPTDF g{t) as function of time t for 
D/Vo = 1.9 and il = 0.2, 0.3, 0.5, 1.2. 
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FIG. 5. Peak areas as a function of Q and MFPT 
< t > as a function of Q. 



In a similar manner, for a fixed (D/ \ Vq \— 1.9), the normalised FPTDF g{t) are plotted for different 
signal frequencies in Fig. 4. It also shows peak structure over the background. As frequency increases the 



peaks are more dominant and for a given frequency peak heights fall exponentially with time. Peaks are 
separated by the period Tq of the signal. After subtracting the exponential background, the peak area of the 
first peak (which is the most dominant) is plotted as a function of frequency of the signal in Fig. 5. It also 
exihibits the non-monotonous behaviour. The curve shows that the synchronization (defining the strength 
of the synchronization as the area of the peak) is maximum at specific value of the signal frequency. 

The non-monotonous behaviour of the peak strength as a function of noise strength or frequency clearly 
indicates that the cooperation between noise and signal becomes maximum at specific value of noise strength 
for given frequency or at specific value of signal frequency for given noise strength. This is the resonant 
behaviour in accordance with our usual convention. Thus this specific value of noise strength and signal 
frequency are denoted as Dres , ^res respectively. 

As MFPT < f > is a very important time scale in this noise induced transition, we also plot < t > after 
evaluating them from Eq.(13) as a function of {D/ \ Vq |) in Fig. 3 and as a function of signal frequency fi 
in Fig. 5. For the sake of comparison MFPT when there is no signal {A — 0), is also plotted as a function 
of (D/ I Vo I) in Fig. 3. The curves show a monotonous fall as noise strength or frequency increases except a 
little curb at low frequency. 
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FIG. 6. MFPT at resonance 
< t >res as a function of period of 
the signal Tq with signal {A — 0.2) 
[line with filled squares] and with- 
out signal (A = 0) [line with filled 
circles! . 
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FIG. 7. Noise strength at reso- 
nance Dres as a function of fre- 
quency Q 
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FIG. 8. Eigenvalues Ai and 
A2 as a function of (| Vo \ /D) 



From the resonance values {Dres/ I K) I) (for the first dominant peak) for different frequencies 17 = 
0.2, 0.4, 0.6, 0.8, 1.O, we evaluate MFPT < t >res from Fig. 3 and plot as a function of period Tq of the signal 
in Fig. 6. The points fit a straight line :< t >= .686T0 -I- 2.2. In order to have a comparison with MFPT 
for these resonant values {Dres/ | Vq |) for non-biased process {A = 0), these values are also plotted as a 
function of the period Tq of the signal in the same Fig. 6. They also fit straightline :< t >~ .384ro + 2.985. 
This figure shows that as period of the signal increases the resonant time scale increases exactly in a linear 
fashion. 

We have seen that synchronization is maximum for specific noise strengths Dres for different signal fre- 
quencies fi. With these resonant values, Dres is plotted as a function of signal frequency fi in Fig. 7. The 
curve is fitted to a straight line [Dres = .42517 + .21], which shows that as frequency increases the noise 
strength for which resonance occurs also should increase linearly. 

In conclusion, we analyse the stochastic synchronization in two-state spatially extended system under the 
influence of periodic field analytically. The analytic expression of the conditional probability could very well 
exihibit the resonance behaviour and synchronization between noise and periodic signal. The expression 
can be further approximated with only the first term in the summation by noting that the second and 
consequently second onward eigenvalues (as the eigenvalues are ordered) differ from it by a large amount. 
The eigenvalues Ai, A2 are plotted as a function of (| Vq | /D) in Fig. 8. The figure clearly shows that for a 
very large variation of the noise strength this approximation is fairly accurate. 

In asymptotic expansion as given in Eq.(3) we consider only the first term. From Eq.(5b) one can estimate 
the magnitude of /i by expanding tpr{x)fi in the complete set of cigcnfunctions {V'n(2;)} at each time t. 
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It can be shown that /i would be a periodic function and independent of position x. In fact, /{ (t) 

1 

[l + (A,/0)2]l/2 



^Im[e^^* / (Xr+ifl)] and consequently | /} (t) |< ri+fA Jnvn/2 < 1- Thus for small (A/D) the approximation 



(7) is fairly good. 

The process (1) has in built scaling properties. When we change time as t -^ Dt, \ Vq \ would change to 
{\Vo\ /D),n-> (n/D), A -^ [A/D) respectively. 

The MFPT of the process at resonance is shown to increase linearly with the period of the signal. The 
resonant noise strength is also shown to increase linearly with frequency of the signal. This fact might help 
the experimentalists in studying the resonance behaviour in complex systems where tuning the frequency is 
a convenient task than varying the internal noise or temperature. 
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